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The dynamical response of Coulomb-interacting particles in nano-clusters are analyzed at different tempera¬ 
tures characterizing their solid- and liquid-like behavior. Depending on the trap-symmetry, both the spatial and 
temporal correlations undergo slow, stretched exponential relaxations at long times, arising from spatially corre¬ 
lated motion in string-like paths. Our results indicate that the distinction between the ‘solid’ and ‘liquid’ is soft: 
While particles in a ‘solid’ flow producing dynamic heterogeneities, motion in ‘liquid’ yields unusually long 
tail in the distribution of particle-displacements. A phenomenological model captures much of the subtleties of 
our numerical simulations. 


Study of nano-systems demands refinement of fundamen¬ 
tal concepts, in addition to their technological promises. For 
example, a sharp transition 01 between the solid and liquid 
turns into a smooth crossover in nano-clusters M- While 
the relaxation in simple solids and liquids produce phonon 
and diffusive modes, particles’ motion in confinements show 
anomalous behavior i7l 4lO(l . Experimental evidences from 
dipolar colloids 0 and dusty plasma 101 . point towards 
a rich interplay of long-range interactions and confinements 
yielding subtle dynamics, and has attracted theoretical atten¬ 
tion i3, 13, 3. Then, how good are these ‘solid’ and ‘liq¬ 
uid’ in traps from a dynamical perspective? We address these 
questions focusing on the dynamical properties of Coulomb¬ 
interacting particles in confinements - with and without sym¬ 
metries. 

The ‘solid’ phase of Coulomb-interacting particles in a 
nano-cluster is called a Wigner molecule (WM) because it 
mimics the physics of a Wigner crystal 0 . The WM ‘melts’ 
into a ‘liquid’ in circular yj and irregular 0] confinements, 
termed circular and irregular Wigner molecule (CWM and 
IWM). A comparison between the two elucidates the role of 
‘disorder’, coming from the irregularities in the trap, on the 
crossover. 

Analysis of melting using spatio-temporal correlations 
yields crucial dynamical information II1611 . The interplay of 
interactions and Brownian motion (consistent with the trap- 
geometry) often leads to frustrated motion and dynamic het¬ 
erogeneities (DH) 10, i.e. multiple time-scale relaxations. 
These are ubiquitous near the transitions in colloids |18|, 
gels 0, glasses 1:23. 0 . and bio-molecule suspensions 1231 . 

We consider simple models for IWM and CWM and an¬ 
alyze the equilibrium motion of the constituents using iso¬ 
kinetic molecular dynamics (MD) simulations 112311 . Our main 
results can be summarized as follows: 

(1) The motion in traps display a variety of non-Gaussian be¬ 
havior of the distribution of particle displacements, includ¬ 
ing its stretched-exponential spatial decay. Such an ultra-slow 
relaxation is consistent with preliminary experimental signa¬ 
tures [;24]. 

(2) The dynamics at long t is intriguing: At low T the solid 
‘flows’ producing DH, whose nature depends on the confine¬ 
ment. Though particle displacements are random at large T, 
their distribution defies central limit theorem 125]. 


(3) Our results demonstrate how the innate dynamics sets 
apart the manifestation of Wigner physics II1511 based on the 
trap symmetries. 


We model a WM using the Hamiltonian: H = <f> G — 

O' I ~~ 1 +E i Konf (g), where the first term represents Coulomb 
repulsion between particles of charge q and <f> = q 2 /e, and 
e is the dielectric constant. Here, I4onf is the confinement 
potential. We compare the results from two confinements: 
(a) Irregula r = a{x A /b + by 4 — 2A x 2 y 2 + j(x — 

y)xyrj 0 0, and (b) Circular V^ T nf (r) = ar 2 , where 
a = 1711x^/2. 


We rescale the length (r — 
cj) 2 / 3 a 1 ^ 3 E), that transform 

H Cr = Ei A + E 

by t 


ji<j 


1/,3 a _1 / 3 r) and energy (E —> 
3(l the CWM Hamiltonian to: 
. The time-scale is set 


h <f~ 2 l 3 a~ x l 3 t [27]. The overall factor a and a 
are chosen to keep mean density in the two traps same for 
a fixed total number of particles, N = 150. The irregularity 
parameters b = n/4, A £ [0.565,0.635] (controls chaotic- 
ity), and 7 £ [0.10,0.20] (breaks reflection symmetry) 1291 
are adjusted to access the standard features of disordered sys¬ 
tems. Distances are measured in the unit of ro, the mean inter¬ 
particle spacing between the neighbors. The statistics of the 
results from V^ nf were improved by ‘disorder averaging ’0 
over five configurations, each identified by a specific {A, 7}. 
For V^ f , averaging is done over five independent MD simu¬ 
lations. 


A thermal crossover from a ‘solid’ to ‘liquid’ in WM as¬ 
serts Si a ‘solid’ at T = 0.006 and ‘liquid’ at T = 
0.03 for our chosen parameters. The spatio-temporal evolu¬ 
tion of the system is studied from the Van Hove correlation 
function (VHCF) ||0 . whose self-part is lEll : G s (r,t) = 
(TV -1 EjLi ^ \ r — |g (t) — g(0)|]). Here (.) denotes average 
over up to 5000 independent time origins. G s (r, t) thus mea¬ 
sures the probability of a particle to move (on average) a dis¬ 
tance r in time t. Focusing on IWM at T = 0.006 in Fig. 
1(a), the dynamics can be broadly classified into three tem¬ 
poral regimes 0: (i) The particles move ballistically for 
t < 1, and the resulting G s (r, t) is Gaussian. They start 
feeling each other only for t > 1 implying that our unit 
time is the time scale of interaction, (ii) Particles get ‘ar¬ 
rested’ in the cage formed by neighbors in the intermediate 
time (100 <t< 1000), leading to a slowly decaying G s [r, t) 
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FIG. 1: The decay of G a (r,t), normalized by G a (r = 0, t) for 
clarity, for several t in the ‘solid’ (T = 0.006): (a) IWM, and (d) 
CWM. G s (r) is Gaussian up to t ~ 1 and displays non-Gaussian 
trends for t > 10, including multi-peaks at large times. Bunching- 
up of the traces of G s (r) for 100 < t < 1000 is shown in the insets 
demonstrating caging (weaker in CWM). Corresponding Af(t) plots 
are shown (b, c: IWM and e, f: CWM) for t = 1000 and t = 9500, 
illustrating the motional coherence. The thick dots represent the ini¬ 
tial position and the connecting lines signify the displacements of the 
same particle at t. Red, green, blue, and purple dots represent par¬ 
ticles with increasing |Ar(f)|.Mobile particles reside in the central 
region for V^ mt , leaving nearly immobile ones in the narrow limbs 
near corners. In contrast, the mobile particles in CWM circle around 
the periphery, while the slow ones flock towards the origin. 


whose traces bunch up, see the inset of Fig. 1(a). (iii) The dy¬ 
namics for t > 1000 develop multi-peak structure in G s (r, t). 

In order to develop a deeper understanding, we next ana¬ 
lyze the particle displacements, Af)(i) = {ry(i) — fy(0)}. 
Most particles only rattles around their equilibrium positions 
for t <C 1000. Motion that begins by breaking a cage causes 
avalanches of such events correlated in space, which leads to 
a string-like path of delocalization B81 I3411 . Fig. 1(b) presents 
Ary (t = 1000) for an IWM at T = 0.006, showing such 
an alignment of neighboring particles leading to a slow relax¬ 
ation that yields a long tail in G s (r, £). Once the collective 
motion subsumes large part of the system, particles undergo 
little relative displacements, but an overall drift, producing 




FIG. 2: The evolution of G a (r,t) in the ‘liquid’ (T = 0.03) of an 
(a) IWM and (b) CWM, shows deviations from the Gaussian behav¬ 
ior (except for small t) in both traps though Ar(t) remains random. 
The nature of the evolution of G a (r , t ) is analyzed in (c) and (d) by 
studying the exponent k(t) (see text), k decays continuously from 
2 to nearly 1 at low T. The circular trap yields 2 < k < 1. k 
keeps degrading with t and T in irregular confinements. The inset in 
(d) shows the Gaussian fit of all the three peaks of a representative 
G a (r, t ) in CWM. t-dependence of MSD at different temperatures 
are shown for irregular (e), and circular (f) traps. The dots repre¬ 
sent MSD by integrating r 2 G a (r,t) (see text), and are shown for 
T = 0.006, 0.03. The perfect match between the solid lines and 
corresponding dots persists for all T. 


secondary peak in G s (r,t ). With further increase in t, the 
string-like path degenerates into sub-clusters, each perform¬ 
ing independent collective motion (Fig. 1(c)), and contribute 
to the multi-peak G s (r , t) of Fig. 1(a). 

The motion in the CWM at T = 0.006 in Fig. 1(d) is qual¬ 
itatively similar to IWM insofar as the temporal regimes are 
concerned, albeit the caging regime is narrow (10 <t< 100) 
and untenable, causing weaker agglomeration of G s (r, t) (in¬ 
set, Fig. 1(d)). However, the differences between the two traps 
stand out in the nature of coherent motion. While the circu¬ 
lar symmetry forces such motion along its azimuth (Fig. l(e, 
f)), the collective dynamics in IWM is guided by the irreg¬ 
ularity producing nonequivalent regions (Fig. l(b, c)). Un¬ 
like the sub-clusters with independent motions in IWM, a sin¬ 
gle (azimuthal) coherence grows in the CWM, see Fig. 1(f). 
This results into a rather broad hump of G s (r, t = 9500) in 
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Fig. 1(d). Interestingly, the mobile particles travel by several 
lattice spacings riding on the coherence even in this ‘solid’ 
phase [351] • Signatures of spatially correlated motion of parti¬ 
cles, have been realized in circular plasma crystals El, and 
also in a flocking transition of a driven granular matter in con¬ 
finements El . 

In the ‘liquid’ (T = 0.03), Ar(f) are large and random. 
Yet, G s (r, t) deviates from Gaussian and differences persist 
in the two traps (Fig. 2(a, b)). While the displacements are 
nearly Gaussian in circular traps, it undergoes a progressively 
slower decline for V^ nf . It is not possible to discern sharp 
temporal regimes in such a ‘melted’ state. 

Next, we quantify the large-r decay of G s (r, t) by fitting it 


to ■ 


„-lr K 


(See supplementary material). Fig. 2 (a: IWM, b: 
CWM) presents the time-evolution of the best fit parameter k. 
For low T, k degrades close to unity (both traps) establishing 
a near-exponential decay for intermediate times. With an in¬ 
crease in T, the long-time value of k rises in however, 

broadly remains within 1 < k < 2. In a striking contrast, k 
decreases monotonically for V^ nf , causing a stretched expo¬ 
nential spatial decay. 

Above analysis does not extend beyond t = 400 due to 
the non-monotonic evolution of G s (r,t) at larger t and low 
T. The monotonic fall of k(t ) in the irregular trap, however, 
continues for large T, and hints towards a saturation only for 
t > 8500. Interestingly, k reaches its minimum for T = 0.02, 
a temperature that marks the crossover between the ‘solid’ and 
‘liquid’ Jl. The traces of k(t ) rises for T > 0.02. The system 
resembles a conventional liquid (k m 2 for large I) only at a 
much larger T ~ 0.25. In contrast, a circular confinement 
yields k —>• 2 at large t for T > 0.05 (See supplementary 
material). 

The classification of temporal regimes is also revealed from 
the mean-square displacement (MSD), (A r 2 (t)), shown in 
Fig. 2 (e: IWM, f: CWM). The low-T evolution of MSD 
is surprisingly flat for IWM showing strong caging. The 
corresponding CWM show a weak Brownian behavior. The 
plateau at intermediate t gradually shrinks with T leading to 
a smooth crossover between ballistic and Fickian motion. It 
follows that, (Ar 2 (f)) = f r 2 G s (r,t)d 2 r, and its estimation 
is shown as dots in Fig. 2(e-f), confirming the ‘Fickian yet 
non-Gaussian’ motion El in the WM at large t. MSD sat¬ 
urates when displacements approach the system size. Its oc¬ 
casional jumps in IWM fades away with T. Such jumps are 
not discernible in a CWM even at the lowest T. IWM, thus, 
resembles a supercooled system El more than what CWM 
does. 

A phenomenological description expresses G s (r, t ) as a 
convolution of Gaussian displacements and the distribution of 
diffusivities, P{D) We generalize this to accommo¬ 

date the multi-peak nature of G s (r, t): 


N v 


P(D p ) 


G s (r, t) = / rf A>773FT7T ex P 


p= 1' 


(Air Dpt) 


0 - r P ) 2 

AD p t 


( 1 ) 


where N p is the number of peaks observed in G s (r, t). By 




FIG. 3: The main panels show the comparison of G a (r, t ) obtained 
from direct MD simulations with its evaluation using Eq. (1). The in¬ 
sets show corresponding P(D) and the parameters t. and T appear as 
legends. The nature of G s (r, t) appears largely controlled by P(D), 
from its rich structures due to heterogeneous diffusion, and from its 
long tail. 


fitting the individual peaks (see the inset of Fig. 2(d)), we 
obtain r p and the peak widths. We determine P{D p ) from 
the statistics of particles whose ( r 2 (t )) falls within the peak 
width corresponding to r p Q. A comparison of Eq. (1) with 
the direct simulation is shown in Fig. 3(a-d). The distinct 
r p signals dynamic heterogeneity for t > 1000. The inset 
of Fig. 3(c) shows a stretched-exponential fall [37 1 of P(D) 
for T = 0.03 in V^ nf (r), whereas, its near exponential decay 
characterizes (r). A full P{D) corresponding to a multi¬ 
peak G s (r,t) displays rich structures (insets of Fig. 3(c-f)), 
and the higher peaks contribute near homogeneous weight to 
P(D). 

The multiple relaxations js also probed by the distinct- 
part of VHCF, G d {r,t ) El: G d (r,t) = (2 [N(N - 
i )] -1 Ei# S[r~\f i(t) — g( 0)|]). It represents the proba¬ 
bility of finding at t a different particle at a distance r from a 
location where there was a particle initially. Fig. 4(a) shows 
the lack of positional order at T = 0.006 in the IWM Q 
through depleted Bragg-peaks ll40ll of G d (r,t -f 0), while a 
CWM displays sharper peaks (Fig. 4(b)) due to the azimuthal 
order. 


We substantiate the heterogeneity defining mobility El as 
(\Afi(t)\). A particle is qualified as ‘fast’ (‘slow’) if its mo¬ 
bility lies within the top (bottom) 7% (10%). The long-time 
accumulation of G d (r —> 0) (Fig. 4(a,b)) arises from the fast 
particles occupying the location from where another particle 
moved out, developing spatially correlated motion. In con¬ 
trast, the slow particles trace out Bragg-peaks. Such distinct 
behavior is illustrated as the inset of Fig. 4(b) for V c ^ f . A 
similar picture emerges from the self-part too - while the slow 
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FIG. 4: Upper panels: Gd(r,t) with r for several t in the ‘solid’, 
depicting the temporal relaxation in (a) IWM and (b) CWM. The 
‘correlation-hole’ in Gd(r « 0, t) fills up with time producing strong 
peak at large t. The inset of (a) identifies the contribution of slow 
(brown) and fast (blue) particles in G 3 (r, t). The inset of (b) shows 
that the fast particles contribute to Gd(r —> 0, t), while the slow ones 
trace the frozen ‘solid’. Lower panels: A strong bond-orientation at 
T = 0.006 reflects in {g&{t)) « 1, for all t for both IWM (panel c) 
and IWM (panel d). Its decay in the melted state (T = 0.03) follow 
an exponential and stretched-exponential t-dependences for circular 
and irregular confinements respectively. 


particles produce the small-r Gaussian behavior, the fast parti¬ 
cles contribute to the long tail of G s (r, t), see inset of Fig.4(a). 

The ‘solid’ state of a WM possesses bond-orientational 
order. Each particle has predominantly six neighbors, and 
each of these when connected to the given particle makes 
an angle of 7 r /3 with the line connecting successive neigh¬ 
bors. The corresponding bond-orientational correlation func¬ 
tion: < 7 6 (t) = (exp[z60(f)]) EH is unity for a perfect Wigner 
crystal. Its decay describes how the propensity of the six co¬ 
ordinated neighborhood corrodes in time. Here 9(t ) is the an¬ 
gle fluctuation of a fixed bond. We present ge(t) in Fig. 4(c) 
and (d) for V^ nf and respectively. A bulk 2D isotropic 
liquid shows an exponential fall of ge(t) with t, signaling the 
loss of this order according to the KTHNY theory 114 lh . While 
our results in Fig. 4(d) broadly capture such a trend in V^ f , 
the ‘liquid’ in V^ n{ (Fig. 4(c)) shows a slower, stretched- 
exponential decay: ge(t) ~ exp [— at 13 ] with /3 ~ 0.6. 

In conclusion, we find non-Gaussian distribution of the par¬ 
ticle displacements, correlated motion of the constituents in 
string-like paths, distinct relaxation of slow and fast particles, 
and stretched-exponential decay of both spatial and tempo¬ 
ral correlations in irregular confinement with Coulomb par¬ 
ticles, whereas the circular traps produce different motional 
signatures. Bulk supercooled liquids, which are also dis¬ 
ordered systems with short-range inter-particle interactions, 
share much of these dynamic footprints. Such similarities hint 
towards these findings to constitute the robust dynamical fea¬ 
tures of disordered systems. It will be interesting to explore 


the response functions [42] in different temporal regimes. Our 
results will provide guidance in designing nano-technology 
devices. 
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We begin this supplementary material elaborating the procedure of obtaining k as a function of t, as presented in Fig. 2(c, d) 
of the main paper. 


EXTRACTION OF THE EXPONENT k(t) 


We fit the MD data of G s (r,t) with two functional dependencies: G| mall (r, f) ~ e r for r < r c , and G l f Ige (r,t) ~ 
e~ lr for r > r c for a range of {/. T}, for which the r-dependence of G s (r, t) has a monotonic evolution. For a given t and T, 
we start our fitting procedure by choosing a rather small value of r c , and fit the entire G s (r, t) curve with either of G| mall (r, t) 
or G 1 s arge (r, t) depending on r smaller or larger than the chosen r c . For each region (r smaller or larger than r c ), we get a value 
for x 2 , measuring the goodness of the fit IS il l in that region. Adding its values from the two regions, we get the total x 2 of the 
entire fit for that specific r c . Now, we increase r c in small incremental steps, and calculate total x 2 for each such r c . Finally, the 
optimal values of r c and other parameters, e.g., c,l,k were obtained by considering the r c for which total yf 1 is minimum. The 
data in Fig. 2(c, d) of the main paper are generated by repeating above procedure for the entire range of {t. T). The adjoining 
Fig- HH depicts one such case, where we highlight the above procedure. 

We emphasize that a value of k = 2 implies that the whole G s (r, t) curve has a single Gaussian fall. While it is also possible to 
get k = 2 from a multiple-Gaussian n fall with different width in different ranges of r, our results of k = 2 do not correspond 
to such possibilities. 



FIG. SI: Illustration of the fitting of G s (r, t) curve with two functional dependencies: G™ all (r, t) ~ e T ’ 2 / c for r < r c , and G large (r, t) ~ 

i fc ^ ^ ^ O • • • 

e~ for r > r c . r c denotes the optimal distance for which total x is minimum. 
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EXPONENT k IN THE LIMIT OF LARGE t AND T 


In the main text, we have presented the evolution of k with data limited to t < 400, and T < 0.03. While extrapolation of the 
same analysis is not possible at large t for the lowest temperatures, we here present k(t) up to t ~ 10000 for 0.02 < T < 0.1. 
The two main features to note here are: 

(a) For V^ nf , the traces of k(t) rises with T for T > 0.02, where as initial increase of T for 0.002 to 0.020 lowered them 
progressively (Fig. 2(a) in the main text). Interestingly, T = 0.020 marks the solid to liquid crossover according to static 
properties BS31 I, for which k attains the minimum value (over all T). 

(b)For the monotonic fall of k(t) for all T persists, and shows a tendency toward saturation for large £(> 8500). 

For KoJif, on the other hand, k approaches a value of 2, indicating Gaussian behavior for large t, and T. The data presented 
in Fig. [S2] has results from a single molecular dynamics simulation resulting into large fluctuations, though we made sure that 
they represent the true situation. 
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FIG. S2: Evolution of the exponent k in the limit of large time ( t ) and temperature (T) for the Irregular (panel a) and Circular (panel b) 
confinements. While, for circular confinement, k(t) is mostly bounded between the values 1 and 2, it goes below 1 and shows a tendency 
towards saturation for iregular confinement for T > 0.02. 


ILLUSTRATION OF CAGING EFFECTS FROM DISPLACEMENT PLOTS 

Fig. IS31 a. b) shows the displacements A r(t) = { f(t) — r(0)} of the particles for irregular and circular confinements, re¬ 
spectively, at intermediate time (t = 100) for T = 0.006, illustrating the caging effect. For both confinements, most particles 
(red dots) perform rattling motion around their equilibrium positions. These particles are ‘trapped’ within the cage formed by 
their nearest neighbors, and such an arrest causes G s (r, t) traces to collapse (inset of Fig. 1(a) in the main text) on each other 
for intermediate times. From the alignment of dots and the connecting lines (Fig. I S3 1 a)). we see that the coherent motion sets 
up randomly due to the absence of spatial symmetry for IWM, while such commencement of the coherent motion in CWM 
(Fig. lS31 b)~) is guided by the circular symmetry. 

Fig. lS3f c. d) shows the displacements of the particles for the two confinements in the same time regime for T = 0.030. Here, 
we see that the displacements are random and relatively large in magnitude, for both the confinements. One can expect such 
random motion at this temperature because of the higher thermal energy. 
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FIG. S3: The displacement A f(t) of the particles are shown at T = 0.006 (panel (a), and (b)), and T = 0.030 (panel (c), and (d)) for IWM 
and CWM at t = 100. Red, green, blue, and brown dots represent initial positions of the particles with increasing magnitude of displacement. 
The length of the connecting black lines correspond to the magnitude of the displacements of the particles at the given time t. 


COMPARISON OF G s {r, t) OBTAINED FROM MD SIMULATION, AND FROM EVALUATION OF Eq.( 1) IN THE MAIN TEXT 

In our paper, we claimed (Fig. 2 (c-f)) that G s (r, t) evaluated from Eq. (1) shows broad agreement with its direct calculation 
from the MD simulation. Here, in Fig. lS41 we quantify such ‘broad’ agreement further, for both the confinements. At T = 0.006, 
while we have already seen G s (r,t) evaluated from MD simulation shows good match with its evaluation from Eq. (1) for all 
time regimes for IWM (Fig. 2 (c) in the main paper), CWM shows small deviation at intermediate times (t = 200, in Fig. 2 (d) 
in the main paper). Here, we substantiate this point in Fig. [S4] (a), and (b) for IWM and CWM, respectively. At intermediate 
times (t = 15, 200), we see good agreement for IWM (Fig. |S4](a)), while small deviation persists for CWM (Fig.[S4](b)). At 
T = 0.030, we found that our phenomenological formalism could not capture the large r trend of the G s (r, t) curve for CWM 
(Fig. 2 (f) in the main paper). Here, from Fig.[S4](c), and (d), we see that, while the match between the two methods of evaluating 
G s (r, t) is not so good for intermediate times for both the confinements, match gets better for IWM as time progresses. 
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FIG. S4: Comparison of G s (r, t) from direct MD simulations with its evaluation from Eq. (1) (see main paper) for IWM (Panel (a),(c)) and 
CWM (Panel (b),(d)) at T=0.006 and 0.030, respectively. Thick dots represent G s (r,t) evaluated from Eq. (1), and solid lines are the same 
evaluated from MD simulation. Insets show the distribution of diffusivities P{D) used to evaluate G s (r, t) for the given t. 


DISTINCT PART OF VHCF AT T = 0.030 


We show the distinct part of VHCF, G,/(r, t), at T = 0.030 in Fig. IS 5 f a, b) for the IWM and CWM respectively. As expected 
in a liquid, the large r Bragg-peaks are already washed out at smallest t. Its relaxation with time smoothens out even the weak 
modulations leaving only the strong peak at r = 0. 



FIG. S5: Gd{r, t) with r for several t in the ‘liquid’, depicting the temporal relaxation in an (a) IWM and (b) CWM. The ‘correlation-hole’ in 
Gd(r ~ 0, f) fills up very quickly for both confinements. 
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SLOW AND FAST PARTICLES IN A CWM 


In order to illustrate the complementary role of the fast and slow particles at low T, we identified in the main text (inset of 
Fig. 3(a)) the sections of a single G s (r, t) contributed by them individually for an IWM at intermediate t. Their independent 
roles can also be seen in circular trap and at large t, which we demonstrate here in Fig.lS6l 



FIG. S6: Illustration of the contribution of slow (green) and fast (blue) particles to G s (r,t) for CWM. The contribution of slow and fast 
particles is consistent with the description given in the main text. 
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